Point symmetries of 3D static plasma equilibrium 
systems: comparison and applications 

(Preprint) 

Alexei F. Cheviakov 

Department of Mathematics, University of British Columbia, Vancouver, V6T 1Z2 Canada 

E-mail: alexch@mast.ubc.ca 

February 7, 2008 

Abstract 

Dynamic plasma equilibrium systems, both in isotropic and anisotropic framework, 
possess infinite-dimensional Lie groups of point symmetries, which depend on solution 
topology and lead to construction of infinite families of new physical solutions. 

By performing the complete classification, we show that in the static isotropic case 
no infinite point symmetries arise, whereas static anisotropic plasma equilibria still 
possess a Lie group of symmetries depending on one free function defined on the set 
of magnetic field lines. The finite form of the symmetries is found and used to obtain 
new exact solutions. We demonstrate how anisotropic axially- and helically-symmetric 
equilibria are obtained using conventional Grad-Shafranov and JFKO equations. 

A recently developed multifunctional automated Maple-based software package for 
symmetry and conservation law analysis is presented and used in this work. 

PACS Codes: 05.45.-a , 02.30.Jr, 02.90. +p, 52.30.Cv. 

Keywords: Plasma equilibria; Lie group; Symmetry; Symbolic computations; Exact so- 
lutions; Grad-Shafranov equation. 



1 Introduction. 

The systems of isotropic Magnetohydrodynamics (MHD) and anisotropic Chew-Goldberger- 
Low (CGL) plasma equations, in particular, their equilibrium reductions, are extensively 
used in physical applications, including the controlled thermonuclear fusion research, geo- 
physics and astrophysics (Earth magnetosphere, star formation, solar activity) and labora- 
tory and industrial applications [HI21|3]. 

MHD and CGL systems, as well as their equilibrium versions, are essentially nonlinear 
systems of partial differential equations, depending in general on three spatial coordinates 
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(x,y,z). Knowledge of physically meaningful exact solutions and analytical properties of 
these systems (such as symmetries, conservation laws, stability criteria, etc.) is of great 
importance for understanding the core properties of the underlying physical phenomena, for 
modeling, and for the development of appropriate numerical methods. 

Common ways of finding exact analytical solutions to such systems include: (1) PDE sys- 
tem reduction by a symmetry group (self-similar or invariant solutions), (2) use of intrinsic 
symmetries of PDEs to generate new solutions from known ones, and (3) use of transfor- 
mations from solutions to other equations. The first approach has been successfully applied 
to the static MHD system in cases of axial and helical symmetry (the well-known Grad- 
Shafranov jUE] and JFKO [B] equations), which resulted in several types of exact solutions 
(e.g. [7j, |H]- [IT] ) . A different approach that makes use of equilibrium solution topology 
(general existence of 2D magnetic surfaces) has been suggested in [T^] . 

A symmetry of a system of DEs is any transformation of its solution manifold into itself. 
A symmetry thus transforms any solution to another solution of the same system. Point 
symmetries are obtained by the general algorithmic Lie method of group analysis ^3]. It is 
equally applicable to algebraic and ordinary and partial differential equations. In a similar 
algorithmic manner, contact and Lie-Backlund symmetries (depending on derivatives) are 
found. Closely connected is the potential symmetry method [H] for discovering nonlocal 
symmetries, and the notion of nonlocally- related potential systems and subsystems |15j . 

It is well-known that Lie group analysis method is capable of detecting not only several- 
parameter, but also infinite-dimensional symmetry groups (depending on an arbitrary func- 
tions, see e.g. [To! H7j .) 

In this paper we perform Lie group analysis and the comparison of point symmetries of 
static isotropic (MHD) and anisotropic (CGL) equations. 

Though the Lie procedure is straightforward and applicable to any ODE/PDE system with 
sufficiently smooth coefficients, it requires extensive algebraic manipulation and the solution 
of (often large) overdetermined systems of linear PDEs. For many contemporary models, 
especially those that do possess non-trivial symmetry structure, such analysis presents a sig- 
nificant computational challenge. For example, for static and dynamic MHD and CGL equi- 
librium systems, the determining equations for the Lie point symmetries split into hundreds 
of linear equations. 

We note that the recent derivation of important infinite Lie groups of symmetries of dy- 
namic MHD [T%1 IT^j and CGL [20] equilibrium systems was done by special methods other 
than the Lie algorithm - the group analysis for these systems has net been previously done 
due to computational difficulty. However later it was shown that Lie algorithm indeed can 
detect those symmetries [TBJ[T7]. 

Modern analytical computation methods, based on Grobner bases [21]- [21] and charac- 
teristic sets j2H]-[2E]> efficiently reduce large overdetermined systems of partial differential 
equations, and in many cases facilitate the complete or partial group analysis. (For a review, 
see |2Z|.) 

One of the most powerful software packages for symmetry-related computations is a REDUCE- 
based set of programs by Thomas Wolf [23- It includes CRACK, LIEPDE and CONLAWi 
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routines for PDE system reduction, search for symmetries, conservation laws and adjoint 
symmetries of DE systems. 

In this work, we introduce a recently developed package GeM for Maple and use it for 
the complete Lie group analysis of the static CGL and MHD systems. The package routines 
perform local (Lie, contact, Lie-Backlund) and nonlocal symmetry and conservation law 
analysis of ordinary and partial differential equations without human intervention (Section 
12.31 ) Corresponding determining systems are automatically generated and reduced, and in 
most cases completely solved. The routines employ a special symbolic representation of the 
system under consideration, which leads to the overall significant speedup of computations. 
This factor is of particular importance in the problems of symmetry/conservation law classi- 
fication for systems with arbitrary constitutive function(s). The module routines can handle 
DE systems of any order and number of equations, dependent and independent variables. 
To the best of the author's knowledge, GeM is the most functional and fast end-user-oriented 
package for Maple. 

The complete Lie group analysis and comparison of symmetries of static CGL and MHD 
systems reveals that when the anisotropy factor is constant on magnetic field lines, the static 
CGL system possesses an infinite symmetry group, involving an arbitrary function constant 
on magnetic surfaces, similarly to the dynamical CGL equilibrium system (V 7^ 0) (201- On 
the other hand, unlike its dynamic version [T£l, the static MHD system does not possess 
infinite symmetries, but only affords scalings, translations and rotations (Section 13.11 ) 

The infinite symmetry group of the static CGL system deserves special attention. It de- 
pends on the solution topology, and reveals an alternative PDE system representation that 
directly connects the static CGL system with the static MHD system, establishing a "many- 
to-one, onto" solution correspondence ( Sections I3.1U3.2I ) The latter has important practical 
value. In particular, a direct way of construction of new CGL equilibria from known MHD 
ones follows. For new CGL solutions obtained from a given MHD one, physical conditions 
and explicit stability criteria are preserved, which makes them suitable for exact physi- 
cal modelling. We also demonstrate how axially and helically symmetric anisotropic static 
plasma equilibria can be constructed using the common form of Grad-Shafranov and JFKO 
equations (Section 14.11 ) (We remark that the equation derived in (30] and referred to by 
authors as "The most general form of the nonrelativistic Grad-Shafranov equation describing 
anisotropic pressure effects" is hardly usable, due to its complicated form, and no nontrivial 
solutions to it are known.) 

In Section T4.21 we present an example of the application of the derived infinite symmetry 
group to construction of a family of new exact static anisotropic (CGL) plasma equilibrium 
(an "anisotropic plasma vortex") in 3D space from a known MHD equilibrium (which is at 
the same time a CGL equilibrium). Other types of exact static MHD solutions known in 
literature may be used to produce families of CGL plasma equilibria in the similar manner. 

Further remarks are presented in Section 
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2 General Lie symmetry analysis of static MHD and 
CGL equations 



2.1 The static MHD and CGL equations. 

The static MHD equilibrium system is obtained from the general MHD equations (e.g. pQ) 
in the time-independent steady case V = and has the form 

curl B x B = grad P, div B = 0, (1) 

Here B is the vector of the magnetic field induction; p, plasma density; and P, plasma 
pressure. 

The static time-independent equilibrium of anisotropic plasmas is a similar reduction of 
CGL equations 1201 

B 2 

(1 — r) curl B x B = grad p± + t grad— + B(B • grad r), div B = 0. (2) 
Here r is the anisotropy factor, 



T = ^b^' P ^ = + V± ^ B ^ h3 = X > 2 ' 3 ' ( 3 ) 



and the pressure P is a tensor with two independent parameters p\\,p±- 
The anisotropic equilibrium system (J2J) is closed with an equation of state [2H 1^2"] 

gradr-B = 0, (4) 

which implies that the anisotropy factor r is constant on magnetic field lines. Other possible 
choices include, for instance, empiric relations (see [HI].) (We remark that the relation (j3J) is 
compatible with dynamic equations of state, such as "double-adiabatic" equations suggested 
in the original CGL paper [2j, since time derivatives identically vanish for all equilibria.) 

It is known that in static MHD equilibria the magnetic field B always lies on 2- 
dimensional magnetic surfaces that span the plasma domain. The only possible exception 
is Beltrami-type configurations curlB = aB, a = const. In a compact domain, magnetic 
surfaces are generally tori [52"] . 

In the general anisotropic equilibrium case, 2D magnetic surfaces are not necessarily 
present. However, as shown below, they always exist in static anisotropic plasma config- 
urations ©-©• 

2.2 The Lie procedure. 

For a general system of I partial differential equations of order p, 
E(x,u,u, ...,u) = 0, 

E = {E\ . . . *E l ), x = (x\ . . . , x n ) e X, u = (u\ . . . , u m ) e U, (5) 

( d k u j \ 
?= hn; — tt^-I j = m) eU k , k = l,...,p, 



k \dx il ...dx ik 



the set of all solutions is a manifold Q in (m+n) - dimensional space XxU, which corresponds 
to a manifold Q 1 in the prolonged (jet) space XxU X U\ x • • • x U p of dependent and 
independent variables together with partial derivatives [TB*] . 

In particular, for the systems (JIJ and the independent variables are cartesian 

coordinates x = (xi, x 2 , X3), and the dependent variables u = (£>i, B 2 , B 3 , P) and u = 
(Bi, B 2 , B 3 ,p±, t) respectively (n = 3, m = 4 and m = 5). 

The Lie method of seeking Lie groups of transformations (symmetries) 

(x'Y = f(x,u,e) (z = l,...,n), , g s 
= ^'(x,u,e) (j = l,...,m) 

that map solutions of (0) into solutions consists in finding the Lie algebra of vector fields 

v = e e (x, u) & + £ 77* ( X , u) ^ + £ d 1)fc (x, u, u) ^ 

+-+ E _^ p (x,u,...,u)^. 



81, 



tangent to the solution manifold fi 1 in the jet space. Here e is a group parameter, or a vector 
of group parameters. 

Here ^^ ip are the coordinates of the prolonged tangent vector field corresponding to the 
derivatives u\ x i ; these are fully determined from (£*, r/ fc ) by differential relations. Thus 
tangent vector fields (J7|) are isomorphic to infinitesimal operators 

* = + »>^- (8) 

i k 

Operators (jHJ) are infinitesimal generators for the global Lie transformation group (JHJ), which 
is reconstructed explicitly by solving the initial value problem 

% = e(f,g), f = ^(r,g), (9) 



/*(£ = 0) = x\ g k (e = 0) = u*. 

To find all Lie group generators (jHJ) admissible by the original system (J3J), one needs to 
solve linear determining equations 

VE(X, U, U, U) | E (x,u,u,...,u)=0 = (10) 

1 p 1 V 

on m + n unknown functions £ l , ?7 fc that depend onm + n variables (x, u). 

The overdetermined system on £ l ,r] is obtained from (fTUj) using the fact that <P,?7 fc do 
not depend on derivatives u\. Thus the system JTUj) splits into an < i(mn + 1) linear 
partial differential equations. For static MHD (CGL) equilibria, this leads to a system of 
133 (253) linear PDEs on 7 (8) unknown functions respectively. Methods of solution of such 
overdetermined systems are discussed in the following subsection. 
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Tangent vector field components £\r] k of the symmetry generator (JJJ) may be allowed to 

depend on derivatives u, u . Such symmetries are known as Lie-Backlund symmetries 

1 P 

Unlike Lie symmetries, they act on the solution manifold VL^ in the infinite-dimensional space 
XxUxU\X---xU p x... of all variables and derivatives. Due to this fact, Lie-Backlund 
symmetries can not be integrated to yield a global group representation similar to (JHJ), and 
can not be used to explicitly map solutions of differential equations into new solutions. 

2.3 "GeM" symbolic package for general Lie symmetry computa- 
tion. 

For a system of several PDEs with several dependent and independent variables, overde- 
termined systems resulting from determining equations can be large - consisting of dozens 
or hundreds of linear PDEs. The simplification and integration of such systems, especially 
in the cases of nontrivial symmetry structure (e.g. infinite symmetries), often presents a 
computational challenge. Analytical computation software is traditionally used to produce 
and reduce large linear PDE systems that arise. 

Algorithms for reduction overdetermined systems based on Grobner bases and charac- 
teristic sets techniques have been successfully implemented j^ZJ- However the author had 
difficulty finding an "end-user" package that does not require human intervention in compu- 
tation, has sufficient functionality (i.e. is capable of computing Lie and Lie-Backlund symme- 
tries of arbitrary ODE/PDE systems of reasonable order, performing symmetry classification 
with respect to constitutive functions of the system, doing conservation law analysis, etc.) 
in conceivable time. 

Available general Lie group analysis software packages, such as Desolv by Vu and Carmi- 
nati, partly have desired functionality, including no user intervention in symmetry compu- 
tation and output in the form of a set of symmetry generators (|7j): however, for the systems 
with several spatial and dependent variables, such as static MHD equilibrium system (J2J), 
the analysis can not be performed with usual computation resources in finite time (though 
the Lie symmetry group finite-dimensional). 

A Waterloo Maple - based package GeM (" General Module" ) has been recently developed by 
the author. The package routines are capable of finding all (Lie, contact and Lie-Backlund) 
local symmetries and prescribed classes of conservation laws for any ODE/PDE system 
without significant limitations on DE order and number of variables, and without human 
intervention. 

The routines of the module allow the analysis of ODE/PDE systems containing arbitrary 
functions. Classification and isolation of functions for which additional symmetries / con- 
servation laws occur is automatically accomplished. Such classification problems naturally 
occur in the analysis of DE systems that involve constitutive functions. The symmetry clas- 
sification may lead to the linearization, discovery of new conservation laws and symmetries 
in particular cases (for example, see [TBI llb|.) 

The " GeM" package [2H] employs a special representation of the system under consideration 
and the determining equations: all dependent variables and derivatives are treated as Maple 
symbols, rather than functions or expressions. This significantly speeds up the computation 
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of the tangent vector field components corresponding to derivatives (up to any order) and 
the production of the (overdetermined) system of determining equations on the unknown 
tangent vector field coordinates. The latter is reduced by Rif [33] package routines (which 
is also done faster with derivatives being symbols, rather than expressions.) 

The reduced system can, in most cases, be automatically and completely integrated by 
Maple's internal PDSolve routine. However (in rare cases) the routine is known to miss 
solutions, therefore it is recommended that the solution of the reduced determining system 
(for symmetry or conservation law analysis) is verified by hand in each case. 

We would like to underline that in the study or classification of symmetries and conservation 
laws, the most important and time-consuming task is indeed the simplification (and in the 
case of a classification problem, the case-splitting) of the overdetermined PDE system of 
determining equations. This task is automatically done by GeM routines. The consecutive 
task of solving the reduced system is in most cases much simpler, and usually can be done 
and verified by hand in reasonable time. 

A brief description of symmetry analysis procedure using the GeM package is given in 
Appendix A. 

3 Comparison of point symmetries of static MHD and 
CGL equilibrium systems 

In the current section we present the complete Lie point symmetry classification of static 
MHD and CGL systems (JTJ, ©-(JU), and demonstrate that the static CGL system possesses 
an infinite-dimensional symmetry group (depending on an arbitrary function), whereas the 
static CGL system does not. 

Global transformation groups corresponding to the discovered symmetry generators are 
given. Special attention is devoted to an infinite-dimensional symmetry group of the CGL 
static equilibrium system, that appears to be related to the infinite-dimensional symmetry 
group of dynamic CGL equilibrium system (201 • (This is not the case for the static MHD sys- 
tem - it does not possess infinite-dimensional symmetry groups, though its dynamic version 
does [HI.) 

The infinite symmetries of the static CGL system ©-0 directly relate it to static MHD 
equations, which is important for applications. Namely, it provides a way of construction of 
exact anisotropic plasma equilibria from known static MHD solutions, allow simple formula- 
tion of famous Grad-Shafranov and JFKO equations for anisotropic plasmas, and more (see 
Section H) 

(The "GeM" package described above has been used for computations.) 
3.1 Symmetry classification 

Theorem 1 (Point symmetries of static CGL and MHD equations) (i)The basis of 
Lie algebra of point symmetry generators of the general static CGL equilibrium system 
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consists of the operators 



d „ d 



Xtrans — 7 Ki~ h K4— ; 

f— f OXi dp± 

k=l 



r\ e\ c\ 

Xrot = (cB 2 + dB 3 )— + (-cB 1 - bB 3 )— + {-dB 1 + bB 2 " 



dB 1 v 1 OJ dB 2 x 1 ^&B 3 

9 9 (9 
+(cx 2 + dx 3 )- h (-cxi - bx 3 )- h (-dxi + bx 2 )^—; 

OXi OX2 ox 3 



fc=i ' \fc=i 



corresponding to translations, rotations and scalings (b, c, d, t, s, w, K A are real param- 

eters.) 

(ii) The static anisotropic CGL plasma equilibrium system (HJ) with equation of the state 
(01) affords an additional infinite- dimensional Lie group of symmetries generated by 



4 + 2(1 -^- B2 ^l' <12) 



where F is a sufficiently smooth function of its arguments. 

(Hi) The Lie algebra of Lie point symmetry generators of static MHD equilibrium system 
(CP is spanned by the infinitesimal X trans , X rot , X^ al , X^ al (with p± replaced by P). 



Proof sequence. 

(i) . The system of determining equations for the static CGL equilibrium system (J2J) consists 
of four equations. After setting all coefficients at derivatives to zero, an overdetermined 
system of 253 equations obtained (using the acgen_get_split_sys () routine.) After sim- 
plification with rif simpO routine, the system is reduced to 64 simpler equations (most of 
which are conditions of the form d^ k /dx m = 0, which are readily integrated by hand to 
produce generators (fTTj) . 

(ii) . For the CGL system (j2J with the equation of state (JH), the picture is more complicated. 
Five determining equations split into an overdetermined system of 199 linear equations. It is 
reduced to 49 simpler equations, which show that rj k are independent of the spatial variables, 
whereas £ fc depend on xi,X2,x 3 linearly, providing rotations and scalings. After factoring 
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out the latter, 15 equations are left, from which tangent vector field coordinates of (|12|) are 
found by straightforward integration. 

(iii). For the static MHD system (0), four determining equations are split into an overdeter- 
mined system of 133 equations After simplification with rifsimpO, the system is reduced 
to 49 simpler equations (most of which in the form of a partial derivative equal to zero), 
which are integrated by hand and give rise to operators X tranS} X rot , X^ al , X^ al . 
This completes the proof of the theorem. □ 

Finite form of point transformations. 

Using the reconstruction formula 0, the global transformation groups corresponding to 
symmetry generators (fTTj) are found: 

1. Xtrans'- translations 

Xi = Xi + KiE, i — 1, 2, 3; p' ± = p± + K A a, B' = B, r' = r; 

2. X rot : SO3 group of 3D space rotations, parameterized, for example, by Euler angles <f>, 9, ip: 

x' = BCD*; B'(x') = 5CT>B(x); p' ± (x') = p ± (x), r'(x') = r(x), 



cos0 sin0 \ / 1 \ / cosip sin^ 

B = I — sin cos (j) ; C = I cos 9 sin 9 ; B = I — sin ip cos ^ 
1/ \ -sin^ cos# / V 1 



seal ^ seal'- ^° independent scalings of coordinates and dependent variables 
Xi = txi, B/ = sBi, p' ± = 2sp±, r = r, % = 1, 2, 3; t, s G R. 

4.X^} al : a group of scalings 

x/ = Xi , B/ = B U 2 = 1,2,3; (p ± + ?pj =c(p ± + ^-Y (1-r)' = G(l-r); C G 



Theorem 2 (The infinite-dimensional symmetry group corresponding to Xqo) T7ie 
static anisotropic CGL equilibrium system possesses an infinite-dimensional set of 

point transformations of solutions into solutions given by formulas 

B x = BM(f), n = 1 - (1 - r)M~ 2 (^), (13) 
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B 2 -B 2 

P±i=P± + g"^, (14) 

These transformations form a Lie group generated by the operator X^ Alty) . Here $ = \I/(r) , 
r G M 3 , is a function constant on magnetic field lines (i.e. on magnetic surfaces when they 
exist.) 

Moreover, the static anisotropic CGL equilibrium system (0)-(0) can be written in the form 



curl (Vl - rB) x (y/l - tB) = grad (p ± + r^-^ 



(15) 



div B = 0, B gradr = 0, (16) 
which implies that p± + rB 2 /2 is constant on magnetic field lines: 

B ■ grad ( P± + r^P) = 0, (17) 



The proof of this theorem is given in Appendix B. 

Remark 1. According to the definition (J3J) of the anisotropy factor r, the parallel pressure 
component p\\ transforms as follows: 

P||i=P±i + B? (l-(l-^^)M- 2 (W) 

Remark 2 . The form of the arbitrary function in the symmetry generator uncovers the 
fact that the quantity p± + rB 2 /2 , together with the anisotropy parameter r, is constant on 
magnetic field lines, and leads to an alternative representation (|15 )) -([16 |) of the static CGL 
equilibrium system ©-(JD). 

3.2 Properties of the infinite symmetries ( 1131) - flTlj ) of the static 
CGL system 

As remarked in the Introduction, both MHD and CGL dynamic equilibrium systems (V ^ 0) 
possess infinite-dimensional symmetries. However it is not so in the static case V = 
0: the MHD system only admits scalings, rotations and translations, whereas the 
static anisotropic (CGL) plasma equilibrium system Q-Q still does possess an infinite- 
dimensional set of symmetries. 

The infinite symmetries f)13|) -(|T^ )l . (fTH|) of the static CGL system constitute a subgroup of 
the infinite-dimensional symmetry group G of dynamic CGL equilibrium system found in 
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[2D] , and share many of their properties. Most important properties of transformations and 
new solutions constructed using them are listed below. 



Structure of the arbitrary function. The transformations (|13| ) - (jl4j) . (|18|) depend on 
the topology of the initial solution they are applied to, namely, on the set of magnetic field 
lines \1/ = const. Here \I/ is an arbitrary sufficiently smooth real- valued function constant 
on magnetic field lines of the original static CGL plasma equilibrium {B,p±,p\\}. In the 
following topologies, the domain of the function ^ = ^(r) is evident: 

(a) . The lines of magnetic field B are closed loops or go to infinity. Then the function 
\[>(r) has to be constant on these lines. 

(b) . The magnetic field lines are dense on 2D magnetic surfaces spanning the plasma 
domain T>. Then the function \&(r) has a constant value in on each magnetic surface, and 
is generally a function defined on a cellular complex (a combination of 1-dimensional and 
2-dimensional sets) determined by the topology of the initial solution {B,px,pii}. 

(c) . The magnetic field lines are dense in some 3D domain A. Then the function \I/(r) is 
constant in A. 

Topology of new solutions. Boundary conditions. Physical properties. New solu- 
tions {Bi,p_i_i,p||i} generated by the transformations (JT3J)-(JTHJ) do not change the direction 
of B, and thus retain the magnetic field lines of the original plasma configuration {B,p_i_,p||}. 
Therefore usual boundary conditions for the plasma equilibrium of the type n • B|gx> = (n 
is a normal to the surface of the plasma domain T>) are preserved. 

If the arbitrary function M(\l/) is separated from zero, then the transformed static anisotropic 
solutions retain the boundedness of the original solution; the same is true about the magnetic 
energy B 2 /2. In models with infinite domains, the free function must be chosen so that the 
new anisotropic solution has proper asymptotics at |r| — ► oo. 

Stability of new solutions. No general stability criterion is available for MHD or CGL 
equilibria. However, for particular types of plasma instabilities, explicit criteria have been 
found. In particular, under the assumption of double-adiabatic behaviour of plasma [2] the 
criterium for the fire-hose instability is [33] 

P||-P±>B 2 , (19) 

(or, equivalently, r > 1), and for the mirror instability - 

If the initial solution {B,pj_,p||} is stable with respect to the fire-hose instability (r < 1), 
then according to (fT3"|) 

1 -n = (1 -r)/M 2 {^) > 0, 

and 7~i < 1, hence transformed solutions retain fire-hose stability /instability of the original 
ones. 
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It can be shown that for every initial configuration {B,p_|_,p||}, there exists a nonempty 
range of values which M 2 (\l/) may take so that the mirror instability also does not occur. 
The proof is parallel to that in Section 2 of 20 . 

Group structure. We consider the set Gc of all transformations (jllSjl. (fTlj) with smooth 
M(ty). Each transformation (|14j) is uniquely defined by a pair {a,//(\E')} as follows: 

M(#) = a exp (#(*)), a = ±1. 

The composition of two transformations (a, -f/"(\l/)) and (j3, Kfif)) is equivalent to a com- 
mutative group multiplication 

(a,H)-(fl,K) = (ap,H + K), 

and the inverse (a, if) -1 = (a, —H). The group unity is e = (1,0). Thus the set Gc is an 
abelian group 

G c = © Z 2 . (21) 

It has two connected components; Ay is the additive abelian group of smooth functions in 
M 3 that are constant on magnetic field lines of a given static CGL configuration. 



4 Applications of infinite-dimensional transformations 
of the static CGL system 

4.1 The Grad-Shafranov and JFKO equations for anisotropic plas- 
mas. 

In 1958, Grad and Rubin [3] and Shafranov j3] have independently shown that the system 
(JTJ) of static isotropic MHD equilibrium equations having axial symmetry (independent of 
the polar angle) is equivalent to one scalar equation, called Grad-Shafranov equation: 

Ar~- + 1>„ + I WW = -r 2 P'(# (22) 
r 

The magnetic field B and pressure P have the form 




p = P(v>) 



(r, ip, z) are cylindrical coordinates. Surfaces ip = const are the magnetic surfaces. 

The JFKO equation is another reduction of the system of static MHD equilibrium equations 
(P), which describes helically symmetric plasma equilibrium configurations, i.e. configura- 
tions invariant with respect to the helical transformations 

z — >• z + 7/1, (p — > ip + h, r — > r, (23) 
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where (r, if, z) are cylindrical coordinates. Solutions to the JFKO equation therefore depend 
only on (r, u), where u = z — 'yip. The equation was derived by Johnson, Oberman, Kruskal 
and Frieman [B] in 1958, and may be written in the form 

*-.!*/' r A + I_M>m + 2jm , (# (24) 



r dr \r 2 + 7 2 r y r 2 + 7 2 ( r 2 + -y 2 ) 2 



The helically symmetric magnetic field is 

■t>h = — e r + B x e z + -t> 2 e0, — — — - — , ts 2 — 5— — 7. — , 

(r, z) are cylindrical coordinates. I(ip) and pressure P{ip) are arbitrary functions. As in 
the Grad-Shafranov reduction, the magnetic surfaces here are enumerated by the values of 
the flux function ip = ip(r,u) = const. 

The direct analogy between the system {1} of static isotropic (MHD) plasma equilibrium 
equations and the alternative representation (|15j) - (jl6|) of the static anisotropic (CGL) equi- 
librium system leads to the following theorem. 

Theorem 3 (Axially and helically symmetric static anisotropic plasma equilibria) 

(i). All axially symmetric solutions of the system of anisotropic (CGL) plasma equilibrium 
equations (H^-CT in 3D space are found from solutions of a Grad-Shafranov equation 

y rr _ *!1 + ^ zz + j($) j'($) = -r 2 iV(^), (25) 
r 

where \l/(r, z) is a flux function, Jfi?) and N(*f>) are some functions of^(r, z), and have the 
form 

r = r(^) < 1, 

^ ~ I r r r v r 6 V ' (26) 

rB 2 rB 2 



p± = N(V) - —, p|| = N(V) + 



2 ' -11 2 

Conversely, any solution to the Grad-Shafranov equation Q25\) corresponds to a family of 
axially symmetric anisotropic plasma equilibria determined by \2b}) . 

(ii). All helically symmetric solution of the system of static anisotropic (CGL) plasma 
equilibrium equations (0)-(^P in 3D space correspond to solutions of the JFKO equation 



r dr \r 2 + 7 2 / r 2 + 7 2 [r 2 4- 7 2 ) 2 



where ty(r,u) is a flux function, J(^) and L(ty) are some functions of^>(r,u), and the heli- 
cally ymmetric static CGL plasma equilibrium is given by 

T = r(*) < 1, 

r 1 f\ , Bp+Bp 1 _ 7^)-^ r _ rJQE)+ 7 g r 

B 2 , , rB 2 



p± = L(*) - pi, = L(tf) + 
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Conversely, any solution to the JFKO equation corresponds to a family of helically 

symmetric anisotropic plasma equilibria determined by \2ty) . 

The theorem statement directly follows from the equivalence of the equations ()15|) - (jl61 1) 
and the system (Q). The equation (J16I 2) provides that the anisotropy factor r = r{ty) is 
constant on axially- (helically-) symmetric magnetic surfaces ^ = const. 

Remark. Several types of solutions to Grad-Shafranov and JFKO equations are available in 
literature. Examples of solutions with physical behaviour can be found in Refs. ff\ ITU l3"B]. 



4.2 Example of an exact solution: an anisotropic plasma vortex. 

Below we use the infinite point symmetries (|13|) -(|14j ) . (J 18)) of static CGL equations to obtain 
a vortex-like exact static anisotropic plasma configuration from a known isotropic one. 

In ^U], Bobnev derived a localized vortex-like solution to the static isotropic MHD equi- 
librium system ((TJ) in 3D space. The solution is presented in spherical coordinates (p, 9, 4>) 
and is axially symmetric, i.e. independent of the polar variable <ft. 1 

The solution has the form 

= e p V(p) cos(#) + e e W{p) sin(0) + e^U{p) sin(0), P v = P - p{p) sin 2 (#), (29) 

where 

Bq,P = const, 

U(p)=X n pV(p), p(p)= 1P 2 V(p), W(p) = -V(p)- P V'(p)/2, 
rr f , rj Vo(2A raP ) - V (2X n R) V (2X n R) 

V{ ~ P) =Bo 1 - V (2X n R) ' 1 = Bo l-V (2\ n R) =C ° nSt > (30) 



V (x) 



sm x cos x 



x 3 x 2 



X n , n = 1, 2, ... is one of countable set of solutions of the equation 

(3 - 4:R 2 X 2 n ) sin(2RX n ) - 6RX n cos(2RX n ) = 0, R = const. (31) 

(In particular, RX\ ~ 2.882, RX2 ~ 4.548, RX3 ~ 6.161.) Bobnev's solution satisfies bound- 
ary conditions 

B\ p=R = 0, P\ p=R = P , B\ p=0 = B e z , (32) 

i.e. the region where the magnetic field is nonzero is a sphere of radius R. The magnetic 
surfaces ^ = const inside the sphere are families of inscribed tori of non-circular section, 
separated by separatrices; the number and mutual position of the families depends on the 
choice of value of R, X n . 



1 Though Bobnev's solution was found by an ad hoc method without using Grad-Shafranov equation, it 
indeed corresponds to some solution of the Grad-Shafranov equation (|22|) . However, it is not possible to 
write the Grad-Shafranov flux function tjj for this solution in the closed form, due to the use of spherical 
coordinates. 
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We use the the point symmetries (fTH |) .()14 p .()18 p . starting with the initial solution (B, r, p±) = 
(B„,0,P„) given by (j2HJ), ((201), corresponding to R = 1,A 3 « 6.161. For this MHD equilib- 
rium, the profile of levels of constant magnetic energy = const and pressure P v = const 
are shown on Fig. 1 (left and right, respectively.) The set of magnetic surfaces \P = P/P m ax 
coincides with levels of constant pressure P v = const. 

As an example, in transformation formulas ([13 |) .([14 p .([18j ) . we choose the arbitrary function 
M(\P) = 1 + ^sin(^), which is separated from zero (since < 1) and constant on the 
magnetic surfaces of the original static MHD configuration. 

Fig. 2 shows the surfaces of constant level of anisotropic pressure components P, p± , p» , 
and their profile along the radius of the vortex in the direction perpendicular to z. 

Fig. 3 illustrates the profiles of magnetic energy densities B^/2 and B^/2 for isotropic and 
anisotropic plasma vortex (left), and the profiles of anisotropic pressure components p± ,p\\, 
in the radial direction (perpendicular to the axis of symmetry z.) 



Magnetic energy levels (Isotropic case) Isotropic pressure P 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: Constant levels of magnetic energy density (left) and pressure (right) in the static 
isotropic MHD equilibrium (Bobnev's solution (j29!30p ) for R — 1, A=A 3 « 6.161. z (upward) 
is the symmetry axis, x (horizontal) is the radial axis. 

Similarly to the isotropic vortex, the anisotropic solution has the boundary condition 
B| p= # = 0, and B| p=0 = B e z at the origin. The magnetic field is chosen to be zero 
outside the plasma domain. 

As seen from Fig. El (right) and expressions (|T4* jl .(|T8 |) . plasma pressure components p\\i,p±i 
of the anisotropic vortex are different within the plasma domain r < 1, but equal on the 
boundary. Thus the presented solution describes a localized anisotropic plasma formation 
confined by external gas pressure. (We note that both known vortex-like isotropic solutions, 
by Bobnev jTUj and by Kaiser and Lortz [32] have the same feature: the plasma ball is 
confined by outer gas pressure.) 
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Anisotropic pressure P perpendicular Anisotropic pressure P parallel 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2: Surfaces of constant level of plasma pressure components P, p± ,p\\ in an 
anisotropic vortex. 




Figure 3: Left: Comparison of energy densities and for isotropic and anisotropic 
plasma vortex in the radial direction. Right: Comparison of plasma pressure components 
Pi Pi- i V\\ i n an anisotropic vortex in the radial direction. 



16 



5 Further remarks 



1. The Lie symmetry-based approach remains to be one of the most generally applicable 
and successful coordinate-independent techniques for generating exact solutions and reduced 
(invariant) versions of nonlinear partial differential equations for obtaining self-similar so- 
lutions. For an ODE, knowledge of Lie symmetries always leads to the reduction of order 
of the equation under consideration, by the number up to the dimension of the symmetry 
group ^3]. The global form of Lie point symmetries is used to generate several-parameter 
(in some important cases, infinite-parameter) families of exact solutions from known ones. 

In this paper we have presented an example when the complete straightforward Lie analy- 
sis reveals an infinite-dimensional symmetry group for static anisotropic (Chew-Goldberger- 
Low) Plasma Equilibrium equations. The existence of this symmetry Lie group leads to a 
dramatic extension of the set of known solutions, bringing in physically relevant anisotropic 
static plasma configurations related to solutions of static MHD equations. In particular, 
one may construct new static CGL equilibria in axial and helical symmetry using corre- 
sponding Grad-Shafranov and JFKO equations ( Section 14. lj) . or from known static MHD 
configurations (as illustrated by the example in Section 

2. The potential symmetry theory ^3] and the algorithmic framework for nonlocally-related 
potential systems and subsystems ^3] has been demonstrated to be useful for calculating 
new nonlocal symmetries and new nonlocal conservation laws for a given system of PDEs, 
especially in the case of two independent variables. However in a PDE system with n > 3 
independent variables, such as MHD or CGL systems discussed in this paper, the corre- 
sponding potential system is underdetermined, and requires suitable gauge constraints (in 
the form of additional equations on the potential variables) to be imposed in order to find 
nonlocal symmetries j3H]- The question of choosing gauge constraints that lead to potential 
symmetries is yet to be answered. Future work will concentrate on this problem, in partic- 
ular, for the system of static MHD equations, which can be totally written as a set of four 
conservation laws [3l)| . 

3. Systems of static MHD ((H) and CGL ©-© equations are examples of systems where 
the use of appropriate analytical computation software is helpful for the full point symmetry 
analysis and leads to finding new symmetry structure (Section |3]) In the same manner, 
infinite symmetries for dynamic MHD [TSj and CGL [2U] equilibrium equations can be found 
within Lie group analysis framework ITS] ITT] , but have been overlooked due to incomplete 
analysis in the preceding literature (e.g. |37|). 

The routines of GeM package for Maple presented and used in this work (and in parts in 
CUE!) perform local (Lie, contact, Lie-Backlund) and nonlocal symmetry and conservation 
law analysis of ordinary and partial differential equations automatically, i.e. without human 
intervention. Corresponding determining systems are automatically generated and reduced, 
and in most cases completely solved. The package routines allow the solution of important 
symmetry/conservation law classification problems for systems with arbitrary constitutive 
functions. To the best of the author's knowledge, GeM is the most complete and fast end- 
user-oriented package for Maple. 

Further development of the GeM package that will include routines for explicit reconstruction 
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of conservation law fluxes and densities and automatic construction of trees of nonlocally- 
related potential systems and subsystems ^5] is one of the directions of current work. 

4. Further research, jointly with G. Bluman and other collaborators, will include the com- 
plete symmetry and conservation law analysis of dynamic CGL and MHD systems, as well as 
the construction of trees of nonlocally-related potential systems and subsystems and search 
for nonlocal symmetries and conservation laws of Nonlinear Telegraph (NLT) equations and 
equations of gas/fluid dynamics (which is partly complete). 
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A "GeM" Symmetry computation module: a brief rou- 
tine description 

Here we outline the program for group analysis of static MHD equations (0) using the GeM 
package. 

First, dependent and independent variables, free functions and free constants are declared. 
For example, for the CGL static equilibrium system (J2J), the statement is 

acgen_decl_vars ( [x , y , z] , [Bl (x , y , z) , B2 (x , y , z) , B3 (x , y , z) , P (x , y , z) ] , [] , [] , 0) ; 

Here Bl, B2, B3 are three components of the magnetic field vector B; the empty bracket 
pairs [ ] indicate that no arbitrary functions and constants are present; stands for the 
highest order of derivative on which arbitrary functions depend. 

Then equations are declared by calling the routine 
acgen_decl_eqs ( 

[ equ[l], equ[2], equ[3] , equ[4] ], 

[diff (Bl(x,y,z) ,x), diff(P(x,y,z),x), diff (P(x,y,z) ,y) , diff (P(x,y,z) ,z)] 
); 

where equ[l] , ...,equ[4] are properly prepared Maple expressions for static MHD equations 
©• 

Within this routine, all derivatives are defined as symbols: dB\jdx = Blx, etc. For 
example, the Maple expression for the equation div B = 

-^Bl (x, y, z) + ^-B2 (x, y, z) + ^-B3 (x, y,z) = 
ox ay az 

is represented as 

Blx + B2y + B3z = 0. 

After this conversion, the unknown tangent vector field components £ l ,r] k are not com- 
posite functions but functions of all "independent variables", which include dependent and 
independent variables of the system under consideration and all partial derivatives. 

The variables on which tangent vector field components if 1 depend are specified by using 

acgen_def_vf ( [Bl(x,y,z) ,B2(x,y,z) ,B3(x,y,z) ,T(x,y,z) ,P(x,y,z)] ,0,0) ; 

The function names of resulting tangent vector field components are stored in a set V_F_C00RDS_. 

The overdetermined system of determining equations system variables on which depend 
are determined by using 

acgen_get_split_sys() ; 
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(which generates 133 equations) and reduced by 



rss : =rif simp (sys_12 , convert (V_F_C00RDS_ , list) ) ; 

to obtain 49 simple equations found in the variable rss [Solved] . 

The final answer for tangent vector field components in this case is easily obtained by hand, 
or Maple built-in integrator 

expand (pdsolve (rss [Solved] ) ) ; 

and provides the results listed in the formulation of Theorem ^ 

We remark again that the built-in Maple routine "pdsolve" for PDE integration sometimes 
produces incomplete answers and therefore is recommended only for preliminary analysis; 
the completeness of its output is to be verified by hand in every case. 

B Proof of Theorem 2. 

Proof of Theorem 2. 

First we prove the formula (|15j). Using the vector calculus identity 

curl(as) = a curl(s) + (grad a) x s 
and the formilas B • gradr = and (0), one gets 



curl (VI - rB) x (VI - rB) = (1 - r) curlB x B - - w (grad T x B) x B 



B 2 B 2 / B 2 

= grad p± + t grad— + — gradr = grad I p± + r— 

We thus observe that the quantity p±+r^- is constant on magnetic field lines: grad (p± + r^- 
B. The same is true for r. Hence without loss of generality 



We may also check that the quantity p± + is invariant under the action of X^: 



where \& is a function constant on magnetic field lines. 

B- 

2 



d ( B 2 V, d . , B 2 d .. d B' 

de [ p± + r T U = dS p ^° + U + T dS— le =° = 0; 



hence ^ may be chosen so that it is also invariant under the action of X^. 
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This allows the integration of Lie equations (JDJ), and the global group of transformations is 
B kA = B k exp (aF (*)) , (1 - n ) = (1 - r) exp (-2eF (*)) , 

P±i = P± + ^ (1 - exp (-2eF (*))) , 

In (|13jl . (|14)l . M is a sufficiently smooth function separated from zero. It is connected with 
F as follows. The case M > is obtained by assigning 

Af(¥) :=exp 

where F is an arbitrary smooth real-valued function; the case M < - the same way by 
taking 

F($) = - + G($), 

where G is an arbitrary smooth real- valued function. 
This completes the proof of Theorem 2.1 
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